BROWNIAN MOTIONS ON METRIC GRAPHS II 
CONSTRUCTION OF BROWNIAN MOTIONS ON SINGLE VERTEX GRAPHS 

VADIM KOSTRYKIN, JURGEN POTTHOFF, AND ROBERT SCHRADER 

Ab STRACT. Pathwise constructions of Brownian motions which satisfy all possible 
boundary conditions at the vertex of single vertex graphs are given. 



1. Introduction and Preliminaries 

This is the second in a series of three articles about the construction and the basic 
properties of Brownian motions on metric graphs. In the first of these articles [21], 
Brownian motions on metric graphs have been defined, their Feller property has been 
shown, and their generators have been determined, i.e., the analogue of Feller's theo- 
rem for metric graphs has been proved. In the present article, we construct all possible 
Brownian motions on single vertex graphs (cf. below). In the third article of this se- 
ries [22], this construction will be extended to general metric graphs. In a companion 
article [20], which serves more as a background for this series and which can also be 
read as an introduction to the topic, we revisit the classical cases [7-9, 14, 15, 18] of 
Brownian motions on bounded intervals and on the semi-line R+. 

For a general introduction to the subject of this article we refer the interested reader 
to [21], henceforth quoted as "article I". We shall refer to equations, definitions, the- 
orems etc. from article I by placing an "I" in front. For example, "formula (1.2.4)" 
refers to formula (2.4) in article I, while "definition 1.3.1" points to definition 3.1 in 
article I. Unless otherwise mentioned, we continue to use the notation and the con- 
ventions set up in that article. 

In the present article we solely consider a metric graph Q consisting of a single 
vertex v, and n,n< +oo, external edges labeled 1^, k e {1,2, . . . , n}, each of which 
is metrically isomorphic to the interval [0, -|-oo), the vertex corresponding to 0. It will 
be convenient to denote the local coordinate of a point in Q° = Q \ {v} by {k, x), 
k e {1,2, . . . , n}, X G (0, +oo), instead of (Ik, x) as we did in paper I. 

The main ideas for the construction of Brownian motions with boundary conditions 
at the vertex compatible with Feller's theorem 5.3 are those which can be found in the 
work by Ito and McKean [14, 15] for the case of the semi-line R_|_, cf. also [18, Chap- 
ter 6]: The reflecting Brownian motion in the case of M_|_ is replaced by a Walsh 
process [31] (cf. also, e.g., [2]) on the single vertex graph, and then the kilhng and 



Date: December 6, 2010. 

2010 Mathematics Subject Classification. 60J65,60J45,60H99,58J65,35K05,05C99. 
Key words and phrases. Brownian motion. Feller Brownian motion, metric graphs. 



2 



V. KOSTRYKIN, J. POTTHOFF, AND R. SCHRADER 



slowing down of this process on the scale of its local time at the vertex is used to 
construct processes implementing the various forms of the Wentzell boundary condi- 
tion. On the other hand, we provide a number of arguments which — at least on a 
technical level — are rather different from those found in the standard literature. For 
example, whenever possible, we use arguments based on Dynkin's formula to derive 
the domain of the generator (i.e., the boundary conditions). This approach appears to 
be much simpler and more intuitive than the one with standard arguments [14, 15, 18] 
for the semi-line R_|_, which is based on the rather tricky calculation of heat kernels 
with the help of Levy's theorem. Moreover, for the case of kilhng, instead of using 
the standard first passage time formula for the hitting time of the vertex we use a first 
passage time formula for the Ufe time of the process. In the opinion of the authors this 
leads to much simpler computations of the transition kernels than those in [14, 15, 18] 
for a Brownian motion on M_|_. 

The article is organized as follows. In several subsections of the present section we 
set up some notation and discuss some preparatory results. In section 2 we recall the 
construction of a Walsh process on a metric graph. In section 3 we construct a Walsh 
process on the single vertex graph with an elastic boundary condition at the vertex, 
while in section 4 we construct a Walsh process with a sticky boundary condition at 
the vertex. The most general Brownian motion on the single vertex graph is obtained 
in section 5 by combining these two constructions. In all three cases we also derive 
explicit expressions of the analogues of the quantum mechanical scattering matrix on 
single vertex graphs. 

1.1. Feller's Theorem and Boundary Conditions. For ease of later reference, in 
this subsection we state Feller's theorem, theorem 1.5.3, as it reads for a single vertex 
graph with n edges. 

Co{Q) denotes the space of continuous functions on Q which vanish at infinity, 
and Cq{G) the subspace of Co{Q) consisting of those functions / G Co{Q) which 
are twice continuously differentiable on Q°, such that /' vanishes at infinity and /" 
belongs to Co{Q). Moreover, for a function / which is continuous on Q°, and such 
that for k G {1,2,..., n}, /(^) has a limit when ^ e Q° approaches the vertex v 
along any edge Ij-, we set 

f{vk)= lim m- 

It is not hard to see that for all / G Cq{Q), k £ {1,2, ... , n}, the limits /'(ufc) of the 
derivatives exist. Theorem 1.5.3 states that the generator A (on Co{G)) of a Brownian 
motion on Q is given by 1/2 times the Laplacean with domain 'D{A) C Cq{Q), 
and that there exist non-negative constants a, c, a ^ 1, and a vector b = {bk, k = 
l,...,n) G [0, 1]", with 



(1.1a) 



n 

a + c + '^bk = 1, 
fe=i 
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and 'D{A) consists exactly of those / for which 



n 



(1.1b) afiv) + -f"{v) = Y,hf{vk) 




1.2. Standard Brownian Motion on the Real Line. The construction of Brownian 
motions on a single vertex graph with infinitesimal generator whose domain consists 
of /'s which satisfy the boundary conditions (1.1) is quite similar to the construction 
carried out for the half-line in [14], [15], [18]. This in turn is based on the properties 
of a standard Brownian motion on the real line, cf., e.g., [10-12, 17, 28, 33], and the 
works cited above. For the convenience of the reader, and for later reference, we 
collect the pertinent notions, tools and results here. 

Let {Qx, X € M) denote a family of probability measures on a measurable space 
{W , A'), and let B = (Bt, t G R_|_) denote a standard Brownian motion defined on 
{Q', A') with Qx{Bq = x) = 1, X G M. It will be convenient to assume throughout 
that B has only continuous paths. Whenever it is notationally convenient, we shall 
also write B{f) for Bt, t > 0. Furthermore, we may suppose that there is a shift 
operator 6 : M+ x O ^ Q, such that for all s,t > 0, Bg o 9^ = Bg+t- 

We shall always understand the Brownian family {B,{Qx, x G M)) to be equipped 
with a filtration J' = {J't, t > 0) which is right continuous and complete for the 
family (Qx, x G M). (For example, J' could be chosen as the usual augmentation of 
the natural filtration of B (e.g., [17, Sect. 2.7] or [28, Sect's 1.4, 111.2]).) 

For any ^ c M, we denote by the hitting time of A by B, 



and we note that for all A belonging to the Borel cr-algebra B{M.) of M., is a 
stopping time with respect to J' (e.g., [28, Theorem 111.2.17]). In the case where 
A = {x}, X G M, we also simply write for H^^y We shall also denote these 

stopping times by H^{A) and H^{x), respectively, whenever it is typographically 
more convenient. The following particular cases deserve special attention. Let x G M. 
Then we have (e.g., [15, Sect. 1.7], [17, Sect. 2.6], [28, Sect's n.3, ni.3]) 



(1.2) 



= inf{t > 0, G A}, 



(1.3) 



Qo{H^ G dt) = Qx{H^ G dt) 



g{t,x)dt, 



t > 0, 



where g is the GauB-kernel 
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(1.4) 



t>o, X eR, 



and 



(1.5) Eo^(e-^^i^) = E^le-^^'o) = g-v^Ai^i^ 

Moreover, for a < x < 6 the law of ^^(,} under is well-known (e.g., [15, 
Problem 6, Sect. 1.7]), and its expectation is given by 

(1.6) E^{Hf^^,y) = {x-a){b-x). 
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Denote by = {Lf , t > 0) the local time of B at zero, where we choose the 
normalization as in, e.g., [28] (and which thus differs by a factor 2 from the one used 
in, e.g., [12, 17]): for a; G M, P^^-a.s. 

(1.7) Lf = \iin^\{{s<t,\Bs\<e]), i > 0, 

£4-0 

and here A denotes the Lebesgue measure. Thus, in terms of its a-potential (cf. [3, 
Theorem V.3.13]) we have 

(1.8) uIb{x) = eJ e""* dLf )= ^e-^l^l, a > 0, x G M, 



which provides an efficient way to compare the various normalizations of the local 
time used in the literature. Slightly informally we can write 

(1.9) Lf= fso{Bs)ds, 

Jo 

where 6o is the Dirac distribution concentrated at 0. is adapted to J, and non- 
decreasing. Moreover, for every a; G M, P3;-a.s. the paths of are continuous, and 
is additive in the sense that 

(1.10) Lf+, = Lf + Lf o et, s, t G M+. 

Similarly as above, we shall occasionally take the notational freedom to rewrite Lf 

as L^{t). 

We will need the following well-known result (e.g., [15, Section 2.2, Problem 3]): 
Lemma 1.1. The joint law of\Bt \ and Lt, t > 0, under Qq is given by 

(1.11) Qo{\Bt\edx, Lf edy) =2^^jlLe-^''+y^^/^Uxdy, x,y>0. 

Let = {Kf, r > 0) denote the right continuous pseudo-inverse of L, 

(1.12) Kf = M{t > 0, Lf > r}, r > 0. 

Note that due to the a.s. continuity of L^ we have a.s. L^^ = r. In appendix B of [20] 
we prove the following 

Lemma 1.2. For any r > 

(1.13) Qo{Kf edt) = ^g{t,r)dt, t>Q, 
and 

(1.14) E^{e-^''r)=e-^', A>0 
holds. 

Moreover, we shall make use of the following lemma, which is similar to results in 
Section 6.4 of [17], and which is proved in appendix B of [20]. 

Lemma 1.3. Under Qq, L^ (^H^_^ ^^y^, x > 0, is exponentially distributed with 
mean x. 
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1.3. First Passage Time Formula for Single Vertex Graphs. In this subsection we 
set up some additional notation which will be used throughout this article. Also we 
record a special form of the first passage time formula, equation (1.2.6). 

Let X be a Brownian motion on Q in the sense of definition 1.3.1 defined on a 
family {n,A,F = {Ft, t > 0),(P^, i G G)) of filtered probability spaces. Let 
be the hitting time of the vertex v. Note that for all ^ ^ Q, P^{Hy < +oo) = 1 (cf. 
the discussion in section 1.3). For A > 0, set 

(1.15) exiO = E^{eM-^H,)) = e-^'^(^'''), ^ e Q, 

where E^{-) denotes expectation with respect to P^. The last equality follows from 
formula (1.5). 

Recall that we denote the natural metric on Q by d. We introduce another symmet- 
ric map dy from Q x Q to defined by 

(1.16) d,{^,7])=d{^,v)+d{v,r]), e, r?Gg, 

which is the "distance from ^ to via the vertex v". Observe that if ^, r] e Q do not 
belong to the same edge, then by the definition of d the equality dy{^, rj) = d{^, rj) 
holds true. 

Next we define two heat kernels on Q by 



n 



(1.17) Pit, ^,V) = J2 I'/c (0 9{t, di^, r,)) (r?), 

k=l 
n 

(1.18) p,{t, ^,V) = Y1 (09{t, dvi^, v)) h, iv), 

k=l 

with t > 0, ^, ?7 G t/. 5 is the GauB-kernel defined in equation (1.4). Hence, in local 
coordinates ^ = {k, x), r] = [m, y), x, y > 0, /c, m G {1, 2, . . . , n}, these kernels 
read 

(1.19) p{t,{k,x),{m,y)) = ^e-^-y)'I^H, 



. '"^ 

(1.20) p,{t,{k,x),{m,y)) = ^e-^-+y? I^t 

V 27ri 

The Dirichlet heat kernel p^ on Q is then given by 

(1.21) P^(.t,^,rj)=p{t,^,r])-pyit,^,rj), t > 0, ^, ry G e?. 

It is the transition density of a strong Markov process with state space ^° U {A} 
which on every edge of g° is equivalent to a Brownian motion until the moment of 
reaching the vertex when it is killed, and A denotes a cemetery state adjoined to Q 
as an isolated point. (Recall from section 1.3 of article 1 that this process is not a 
Brownian motion on Q in the sense of definition 1.3.1.) 

The Dirichlet resolvent = (i?^, A > 0) on ^ is defined by 

(1.22) R^fiO = E^{jJ^ e~^'f{Xt) dt), A > 0, ^ G / G B{g). 
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It is easy to see that has the following integral kernel on Q 

(1.23) rf (e,r?) = rxi^,v) - r,,xi^,v), ^.n^Q, 

where for ^,r)eQ, 



" -V2\d{^,v) 

(1.24) r,(e,r?) = 1^,(0 l/,(r/), 

and 



(1.25) 



" -v'2Ad„(C,r;) 

k=i 



ik=i 

In particular, is the Laplace transform of the Dirichlet heat kernel (1.21) at A > 0. 
Recall the first passage time formula (1.2.6): 

where S is any P^-a.s. finite stopping time relative to J". The choice S = Hy gives 
the following result. 

Lemma 1.4. Let X be a Bwwnian motion on Q with resolvent R = {Rx, A > 0). 
Then for all X> 0, ^ eG, f E B{g), 

(1.26) Rxf iO = RxfiO + exiO Rxf{v) 

holds true. 

The following notation will be convenient. For real valued measurable functions 
/, g on Q, with restrictions fk, gk< A; £ {1, 2, ... , n}, to the edges Ik we set 



» n 

if, 9)= / f{09{0dC = J2(fk,9k), 



k=i 

where the integration is with respect to the Lebesgue measure on Q, and 



poo 

{fk,9k)= fk{x)gk{x)dx, 
Jo 



whenever the integrals exist. 

Assume that / G Co{G)- Then for A > 0, Rxf belongs to the domain of the 
generator of X, and therefore to Cq (cf. subsection 1.1). It is straightforward to 
compute the derivative of the right hand side of formula (1.26), and we obtain the 
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Corollary 1.5. For every Brownian motion X onQ with resolvent R = {R\, A > 0), 
and all f £ Co{Q), 

(1.27) {Rxfy{vk) = 2{ex,k,fk)-V2XRxf{v), k E {1,2, . . . ,n}, 
holds true. 

1.4. The Case 6 = 0. The case, where all parameters bk, k = 1, . . . ,n, in equa- 
tion (1.1) vanish, is trivial in the sense that the associated Brownian motion can be 
constructed by a stochastic process living only on the edge where it started, and there- 
fore it is just a classical Brownian motion on M_|_ in the sense of [18, Section 6.1]. 
This case is also discussed briefly in [18], but for the sake of completeness we include 
it here in somewhat more detail than in [18]. 

Consider a standard Brownian motion on M as before, and without loss of general- 
ity assume in addition that the underlying sample space is large enough such that all 
constant paths in M can be realized as paths of the Brownian motion. Construct from 
the Brownian motion a new process by stopping it when it reaches the origin of M, and 
then kill it after an exponential holding time (independent of the Brownian motion) 
with rate /? > 0. We shall only consider starting points x € R+. If /3 = 0, then the 
process is simply a Brownian motion with absorption at the origin. For example, it 
follows from Theorem 10.1 and Theorem 10.2 in [5] that for every /3 > this process 
is a strong Markov process, and obviously it has the path properties which make it a 
Brownian motion on IR+ in the sense of [18, Section 6.1]. Thus, if ^ G ^, ^ G Ik, 
k = 1, . . . ,n, then we just have to map this process with the isomorphisms between 
the edges Ik, k = 1, . . . , n, and the interval [0, -|-oo) into Q to obtain a Brownian 
motion on G with start in ^, such that it is stopped when reaching the vertex, and then 
is killed there after an exponential holding time with rate /3 > 0. 

Let = {Uf, t > 0) denote the semigroup associated with this process. It is 
obvious that for / G Co{g) we get U^fiv) = exp{-l3t)f{v), t > 0. Thus for the 
corresponding resolvent R'^ = {R^, A > 0), and / G Co(^) one finds 

(1.28) >^Rxfiv) - f{v) + mlfi.^) = 0> A > 0. 

Let be the generator of this process, and recall from lemma 1.5.2, that for all 
/ G V{A^), A^fiC) = 1/2 /"(^), ^eg. But then the identity XR^ = A^R^ + id 
implies the following formula 

l{Rlf)"{v) + PRlf{v) = (}. 

For every A > i?^ maps Co{g) onto P(^°). With the choice a = {1 + /S)"^ /3, 
c = (1 + this shows that the process realizes the boundary conditions of equa- 
tions (1.1) with bk = 0, k = 1, . . . ,n. 

Moreover, we can now use equation (1.26) combined with formula (1.28), to obtain 
the following explicit expression for the resolvent with / G Co(^), A > 0: 

(1.29) RlfiO = R^fiO + ^ e-^'^^^'^) f{v), ^ G g. 
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where, as before, Ri^ is the Dirichlet resolvent. 

In order to compute the heat kernel associated with this process on Q, we invert the 
Laplace transforms in equation (1.29). For the first term on the right hand side this is 
trivial, and gives the Dirichlet heat kernel p^, cf. equation (1.21). The second term 
could be handled by a formula which can be found in the tables (e.g., [6, eq. (5.6. 10)]). 
But this formula involves the complementary error function erfc at complex argu- 
ments, and does not yield a very intuitive expression. Instead, we can simply use the 
observation that t H> cxp(— /3t) is the inverse Laplace transform of A (/3 + A)~^. 
Moreover, the well-known formula for the density of the hitting time of the origin by 
aBrownian motion on the real line (e.g., [15, p. 25], [17, p. 96], [28, p. 102]) gives us 
the following expression for the density of the first hitting time of the vertex 

(1.30) P^{H, e ds) = e-'^(«''')'/2^ ds, s > 0. 
Using the well-known Laplace transform (e.g., [6, eq. 4.5.28]) 

(1.31) / e-^^-==e-" /2^ds = e-^", a > 0, A > 0, 

we infer that the inverse Laplace transform of the exponential in (1.29) is given by 
P^{H-t, G ds). Thus we obtain the following heat kernel 

(1.32) 

pO(i,^,dr/) =p^(t,C,r?)d77- e-'5(*-^)P5(if, G ds)) e^idrj), 

= p^it, ^)dr,-( f e-/^(*-) ^ e-'^«'^)V2« ds) e^v), 
with (,r]eQ,t>0, and is the Dirac measure at the vertex v. 

1.5. Killing via the Local Time at the Vertex. We recall from remark 1.3.2, that 
we may and will consider every Brownian motion X onQ with respect to a filtration 
T = {J^t, t>0) which is right continuous and complete relative to (Pj, ^ G ^), and 
such that X is strongly Markovian with respect to T. 

In this subsection we suppose that X is a Brownian motion on the single vertex 
graph Q with infinite hfetime, and such that the vertex is not absorbing. This entails 
(cf. section 1.3) that X leaves the vertex immediately and begins a standard Brownian 
excursion into one of the edges. Therefore we get in this case for the hitting time of 
the vertex Pv{Hy = 0) = 1, i.e., v is regular for {v} in the sense of [3]. Consequently 
X has a local time L = [Lt, i > 0) at the vertex (e.g., [3, Theorem V.3.13]). Without 
loss of generality, we suppose throughout this subsection that L is a perfect continuous 
homogeneous additive functional (PCHAF) of X in the sense of [33, Section III.32]. 
That is, L is a non-decreasing process, which is adapted to T, and such that it is a.s. 
continuous, additive, i.e., Lt+s = Lt + Lg o 6t, and for all ^ ^ G, -P^ (-^o = O) = 1 
holds true. Moreover we may and will assume from now on that X and L are pathwise 
continuous. 



BROWNIAN MOTIONS ON METRIC GRAPHS II 



9 



Killing X exponentially on the scale of L, we can construct a new Brownian mo- 
tion X onQ. We shall do this using the method of [17, 18]. 

Let K = [Ks, s G M+) denote the right continuous pseudo-inverse of L: 

(1.33) Ks = mf{t >0, Lt> s}, s G M+, 

where — as usual — we make the convention that inf = +oo. The continuity of L 
entails that for every s G M_|_, L^^ = s. Clearly, K is increasing, and due to its right 
continuity it is a measurable stochastic process. Fix s G M+. It is straightforward to 
check that for every t G M+, 

(1.34) {Ks <t} = {Lt > s}. 

Because L is adapted, the set on the right hand side belongs to Jj, and since is 
right continuous, equality (1.34) shows that for every s G M+, Kg is a stopping time 
relative to T. We remark that since L only increases when X is at the vertex v, the 
continuity of X implies that for every s G M+ we get X{Ks) = v on {Kg < +oo}. 
On the other hand, we shall argue below that L a.s. increases to +oo, so that we get 
X{Ks) = V a.s. for all s G M+. 

Let /3 > 0. Bring in the additional probability space (R+, i3(M+), P^), where 
is the exponential law with rate /3. Let S denote the associated coordinate random 
variable S{s) = s, s G M+. Define 

A = A(E)B{M.+), 

We extend X, L, K, and S in the canonical way to these enlarged probability spaces, 
but for simplicity keep the same notation for these quantities. 
Set 

(1.35) C/3 = mf{t>0, >5}, 

and observe that since K is measurable we may write = Ks- Thus as above we 
get X{(j3 ) = V. Define the killed process 

(1.36) Xt = 

Since this prescription for killing the process X via the (PCHAF) L is slightly 
different from the method used in [3, 33], we cannot use the results proved there to 
conclude that the subprocess X of X is still a strong Markov process. However, it 
has been proved in [20, Appendix A] that the strong Markov property is preserved 
under this method of killing, i.e., X is a strong Markov process relative to its natural 
filtration (actually relative to a larger filtration, but we will not use this here). Now we 
may employ the arguments in section 2.7 of [17], or in section 111.2 of [28] to conclude 
that X is a strong Markov process with respect to the universal right continuous and 
complete augmentation of its natural filtration. 
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It is clear that X has a.s. right continuous paths which admit left limits, and that its 
paths on [O,^^) are equal to those of X, and thus are continuous on this random time 
interval. Moreover, it is obvious that for every ^ G Q°, we have P^{C0 > H^) = 1. 
Therefore, up to its hitting time of the vertex, X is equivalent to a Brownian motion 
on the edge to which ^ belongs, because so is X. Altogether we have proved that — 
under the hypothesis that L is a PCHAF, which will be argued below in all cases that 
we consider — X is a Brownian motion on ^ in the sense of definition 1.3.1. 

There is a simple, useful relationship between the resolvents R and R of the pro- 
cesses X and X, respectively. Recall our convention that all functions f on Q are 
extended to ^ U {A} by /(A) = 0. 

Lemma 1.6. For all A > 0, / G B{g), ^ e Q, 

(1.37) Rxf iO = RxfiO - exiO {e'^^') Rxf{v) 

holds true, where ex is defined in equation (1.15). 

Proof. For A > 0, / e B{g), ^eQ 

Rxf{i) = E^[j\-^'f{Xt)dt) 

= Rxfi^) - E^{e-^^^ ^"e-^*/(^t+C.) dt). 
By construction, the last expectation value is equal to 

^ i) ^~^J^ e-^'E^(^e-^''^f{Xt+K.))dtds, 

where we used Fubini's theorem. Consider the expectation value under the integrals, 
and recall that for fixed s € M_|_, Ks is an J^-stopping time, while X is strongly 
Markovian relative to J^. Hence we can compute as follows 

E^[e-^^^ f{Xt+K.)) = E^[e-^^^ E^[f{Xt+K. \ Tk^))) 

= E^{e-^^^)E4fiXt)), 
where we used the fact that a.s. X{Ks) = v. So far we have estabhshed 

RxfiO = RxfiO - Rxfiv). 

In order to compute the expectation value on the right hand side, we first remark that 

because L is zero until X hits the vertex for the first time, we find that for given 
s G M+, Kg > Hy, and therefore Kg = Hy + Kg o Ojj^. Hence, and again by the 
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strong Markov property, 

= E^{e-^^^E^{e-^''^°''^\TH^)) 
= E^{e-^^^)E,{e-^^^), 

Integrating the last identity against the exponential law in the variable s, we find with 
formula (1.15) 

and the proof is finished. □ 

Remark 1.7. Formula (1.37) is quite useful, because if the resolvent of X is known, 
then — in view of equation (1.15) — it reduces the calculation of Rx to the computa- 
tion of the Laplace transform of the density of under P^. 

2. The Walsh Process 

The most basic process — which on a single vertex graph plays the same role as a 
reflected Brownian motion on the half line — is the well-known Walsh process, which 
we denote by = {Wt, t > 0). It corresponds to the case where the parameters a 
and c in the boundary condition (1.1) both vanish. This process has been introduced 
by Walsh in [31] as a generalization of the skew Brownian motion discussed in [15, 
Chapter 4.2] to a process in which only moves on rays connected to the origin. 

A pathwise construction of the Walsh process in the present context is as follows. 
Consider the paths of the standard Brownian motion B = {Bt, t > 0) on M, and 
its associated reflected Brownian motion \B\ = {\Bt\, t > 0), where | • | denotes 
absolute value. Let Z = {t > 0, Bt = 0}. Then its complement Z'^ is open, and 
hence it is the pairwise disjoint union of a countable family of excursion intervals Ij = 
{tj,tj-^.i), j G N. Let R = {Rj, j G N) be an independent sequence of identically 
distributed random variables, independent of B, with values in {1, 2, . . . , n} such that 
Rj, j G N, takes the value k G {1, 2, . . . , n} with probability G [0, 1], J2k "^k = 
1. Now define Wt = vift Z, and if t G Ij, and Rj = k set Wt = {k, \Bt\). In 
other words, when starting at ^ G Q°, the process moves as a Brownian motion on the 
edge containing ^ until it hits the vertex at time H^, and then W performs Brownian 
excursions from the vertex v into the edges Ik, k G {1, 2, . . . , n}, whereby the edge 
Ik is selected with probability Wk- 

As for the standard Brownian motion on M (cf. subsection 1.2), we may and will 
assume without loss of generality that W has exclusively continuous paths. 

Walsh has remarked in the epilogue of [31], cf. also [2], that it is not completely 
straightforward to prove that this stochastic process is strongly Markovian. A proof 
of the strong Markov property based on Ito's excursion theory [13] has been given 
in [29, 30]. A construction of this process via its Feller semigroup can be found in [2] 
(cf. also the references quoted there for other approaches). 
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Next we check that the Walsh process has a generator with boundary condition at 
the vertex given by (1.1) with a = c = 0. Let / G V{A'^). At the vertex v Dynkin's 
form for the generator reads 

Ejf{{X(HZ,))) - /(v) 



(2.1, A-m=^ 

where if^^ is the hitting time of the complement of the open ball Bf{v) of radius 
e > around v. 

Lemma 2.1. For the Walsh process Ey{H^^^) = e^. 

Proof. Since by construction W has infinite lifetime, H^^ is the hitting time of the 
set of the n points with local coordinates (k, e), k = 1, . . . , n. Therefore, by the 
independence of the choice of the edge for the values of the excursion, it follows that 
under Py the stopping time H^ ^^ has the same law as the hitting time of the point e > 
of a reflected Brownian motion on M+, starting at 0. Thus the statement of the lemma 
follows from equation (1.6). □ 

From the construction of W we immediately get 



n 



k=l 

with the notation fk{x) = f{k,x), x G M+. Inserting this into equation (2.1) we 
obtain 

n 

^ k=l 

and since f'{vk) exists (cf. section 1.5) it is obvious that this entails the condition 

n 

(2.2) Ywkf{vk) = 0. 

k=l 

For later use we record this result as 

Theorem 2.2, Consider the boundary condition (1.1) with a = c = 0, and h G 
[0, 1]". Let W be a Walsh process as constructed above with the choice = ^fe, 
k G {1, 2, . . . , n}. Then the generator A"' ofW is 1/2 times the second derivative on 
Q with domain consisting of those f G Cq{Q) which satisfy condition (1.1b). 

For the remainder of this section we make the choice a = c = 0, Wk = bk, 
k G {l,2,...,n} in(l.l). 

Next we compute the resolvent of W. Let A > 0, / G Co{Q), and consider first 
^ = V. Without loss of generality, we may assume that W has been constructed 
pathwise from a standard Brownian motion B as described above, and that B is as in 
subsection 1 .2. Then we get 



EvifiWt)) = Y.bmE^{f^{\Bt\)). 



m=l 
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Hence we find for the resolvent of the Walsh process 

(2.3a) Rxf{v)= I r^iv,v)f{r])dv 

Jg 

with resolvent kernel r^{v,r]), rj e Q, given by 

„-^/2\d{v,'r]) 

(2.3b) r^iv,v) = Yl '^bm ' ^ luir^), 

m=l ^ 

and where the integration in (2.3a) is with respect to the Lebesgue measure on Q. 

Now let ^ ^ Q. We use the first passage time formula (1.26) together with formu- 
lae (1.15) and (2.3), and obtain 

Lemma 2.3. The resolvent of the Walsh process on Q is given by 

(2.4a) RlfiO = f rj;{C,v)f{ri)dr,, A > 0, ^ e a, / G B{g), 

JQ 

with 

(2.4b) rl{^, r,) = r,(^, v) + J2 ^^MO ^fem 4^ e,,^(r/), 

k,m=l 

(2 Ac) S^^ = 2Wm - hm. 

where r\ is defined in equation (1.24), and where e\k> ^\,m denote the restrictions of 
e\ (cf. (1.15)) to the edges Ik, Im respectively. 

Remark 2.4. The matrix S"^ = {S'k„^, k,m = 1, . . . ,n) is the scattering matrix 
as defined in quantum mechanics. We briefly recall its construction in the present 
context, for more details the interested reader is referred to [23]. S'^ is obtained from 
the boundary conditions at the vertex v in the following way. Consider a function 
f on Q which is continuously differentiable in ^° = ^ \ {v}, and such that for all 
A; = 1, . . . the Umits 

F;, = f{vk)= lim /'(O 

exist. Define two column vectors F, F' eC^, having the components Fk and Fl^,k = 
1 , . . . , n, respectively. Furthermore, consider boundary conditions of the following 
form 

(2.5) AF + BF' = 0, 

where A and B are complex n x n matrices. The on-shell scattering matrix at energy 
E > Ois defined as 

(2.6) Sa,b{E) = -(^ + i^B)-\A - iVEB), 
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which exists and is unitary, provided the nx2n matrix {A, B) has maximal rank (i.e., 
rank n) and AB'^ is hermitian. These requirements for A and B guarantee that the cor- 
responding Laplacean is a self-adjoint operator on L'^{Q) (with Lebesgue measure). 
Observe that under these conditions the boundary conditions (2.5) are equivalent to 
any boundary conditions of the form CAF + CBF' = where C is invertible. Also 
Sca,cb{E) = Sa,b{E) holds true. For the Walsh process at hand, concrete choices 
for A and B are given by 
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Then (2.5) is the condition that / is actually continuous at the vertex v, i.e., f{vk) = 
f{vm), k, m = 1, . . . ,n, and that (2.2) is valid (with Wk = 6fc, A; G {1, 2, . . . , n}). 
Obviously, {A^, B^) has maximal rank. However, A^{B^)'^ is hermitian if and only 
if all bk are equal (i.e., hk = 1/n, k = 1, . . . ,n). Nevertheless, (2.6) exists also in the 
non-hermitian case, and Sa^,b^{E) = holds for all E > due to the relations 
j^wgw ^ _j^w^ Qwgw ^ jjj addition, the following relations are valid: 

(2.7) 5'" = (5^)"\ 

(2.8) detS"' = (-1)"+^ 

Furthermore, S"" is a contraction, and the associated Laplace operator is m-dissipative 
on (^) since trivially Im(^i3^ ) = 0, cf . Theorem 2.5 in [ 1 9] . When all bk are equal, 
such that ^"'(S"')t = 0, then S"" is an involutive, orthogonal matrix of the form 

(2.9) 5"' = -1 + 2P„. 

Pn is the matrix whose entries are equal to 1/n. Pn is a real orthogonal projection, 
that is Pn = pI = P^ = P^. It is also of rank 1, that is dimRanP^ = 1. The 
relation (2.4b) giving the resolvent in terms of the scattering matrix is actually valid 
in the more general context of arbitrary metric graphs and boundary conditions of the 
form (2.5), see [19,24]. 

It is straightforward to compute the inverse Laplace transform of the right hand 
side of formula (2.4b), and this yields the following result. 
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Lemma 2.5. For t > 0, ^, r] E Q the transition density of the Walsh process on Q is 
given by 

n 

(2.10) p"'(t,^,r/)=p^(t,C,r7)+ 1^.(6 2^m5(i,d,(e,r?)) 1,^(77), 

k.m=\ 
n 

(2.11) =p{t,i,ri)+ Yl h.iOS]:^9{t,d,{^,v))lUv)- 

k,m=l 

pit, ^, 77) is defined in equation (1.17), ^, rj) in equation (1.21), g is the Gaufi- 
kernel (1.4), and is defined in equation (1.16). 

Remark 2.6. Alternatively p'^{t, ^, 77) can also be written as 

P''{t,^,ri)=p{t,^,v) 

^^■^^^ + E l'.(^) fpdHveds)S]^^g{t-s,d{v,v))lUri)- 

Even though this formula appears somewhat more complicated than (2.1 1), it exhibits 
the role of the scattering matrix 5"", that is, it describes more clearly what happens 
when the process hits the vertex. 

3. The Elastic Walsh Process 

In this section we consider the boundary conditions (1.1) with < a < 1 and c = 
0. The corresponding stochastic process, which we will denote by W^, is constructed 
from the Walsh process W of the previous section in the same way as the elastic 
Brownian motion on is constructed from a reflected Brownian motion (cf., e.g., 
[14], [15, Chapter 2.3], [18, Chapter 6.2], [17, Chapter 6.4]). 

In more detail, the construction is as follows. Consider the Walsh process W as 
discussed in the previous section. We may continue to suppose that W has been 
constructed pathwise from a standard Brownian motion B, as it has been described 
there. But then the local time of W at the vertex, denoted by L"^ , is pathwise equal 
to the local time of the Brownian motion at the origin (and we continue to use the 
normahzation determined by (1.8)). It is well-known (e.g., [15, 17, 18,28]) that U" 
has all properties of a PCHAF as formulated in subsection 1.5 for the construction 
of a subprocess by killing W at the vertex. We continue to denote the rate of the 
exponential random variable S used there by ^ > 0. Let be the subprocess so 
obtained. In particular (cf. 1 .5), is a Brownian motion on Q, and in analogy with 
the case of a Brownian motion on the real line we call this stochastic process the 
elastic Walsh process. We write C/3,0 for the lifetime of (i.e., for the random time 
corresponding to C,p in subsection 1.5). 

We proceed to show that the elastic Walsh process has a generator with 
domain V^A^) which satisfies the boundary conditions as claimed. In other words, we 
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claim that there exist a G (0, 1) and 6^ G (0, 1), A; G {1, 2, . . . , n}, with a + bk = 
1, so that for all / G V{G), 

n 

(3.1) af{v) = ^bkf'{vk) 

k=l 

holds. To this end, we calculate A^f{v) in Dynkin's form. We shall use a notation 
similar to the one used in subsection 1.5. Namely, let and Ey denote the probability 
and expectation, respectively, on the probability space extended by (M_|_, i3(M_|_), P^), 
while the corresponding symbols without • are those for the Walsh process without 
kilUng. 

For e > and under let ^ denote the hitting time of the complement B^{vY 
of the open ball i?e(v) of radius e > with center v. Then iJ^ ^ = if^^ A C,^^, 
where as before iJ^g is the hitting time of B^{vy by the Walsh process W. (Note that 
B^{vy contains the cemetery A.) We find 

n 

(3.2) ^,(/(T^^(i?„%))) = Pv{K,e < Cao)- 

k=l 

The probability in the last expression is taken care of by the following lemma. 
Lemma 3.1. For all e, /3 > 0, 

(3-3) PvKe < C0,O) = 

Proof. We may consider the Walsh process W as being pathwise constructed from 
a standard Brownian motion B on the real line as in the previous section, and we 
shall use the notations and conventions from there. Then it is clear that under P„ 
and under P^,, HJ^^ has the same law as the hitting time of the point e in M+ by the 
reflecting Brownian motion \B\ under Qo, that is, as HP of the Brownian motion 
B itself under Qq. Let K'^ denote the right continuous pseudo-inverse of L^. For 
fixed s G M+ we get 

{K^ < H^,} = {L-(if- ) > 

Hence 

P,{K^ <Hi:)=P,{L^{Hl,)>s) 

= Qo{L^{Hf_^^,^)>s). 

In appendix B of [20] is shown with the method in [17, Section 6.4] that under Qq the 
random variable L^{HP_^ ) is exponentially distributed with mean e. So we find 
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We integrate this relation against the exponential law with rate /3 in the variable s, and 
obtain 

/•OO 

Jo 

1 



l + e/3 

We used the fact that due to the continuity of the paths of W we have (/^^ / H^,e- ^ 
We insert formula (3.3) into equation (3.2), and obtain 

= '3S i;^ mm'(K.))) - /(")) 

1 1 " 
lim — ( > Wkfhie) — fiv)] 



k=l 



lim 



1^ ^fe /^/(^) ) ■ 



Obviously Ey{H^,^J < Ey{H^^J = (cf. lemma 2.1). Since the last limit and 
f'{vk), k e {1,2, . . . , n}, exist and are finite, we get as a necessary condition 

n 

(3.4) ^Wkf'{vk)-/3f{v)=0. 

k=l 

Thus we have proved the following theorem. 

Theorem 3.2. Consider the boundary condition (1.1) with a G (0, 1), b G [0, 1]", 
and c = 0. Set 

(3.5) Wk = k = l,...,n, p=—^, 

1 — a 1 — a 

and let be the elastic Walsh process as constructed above with these parameters. 
Then the generator ofW^ is 1/2 times the second derivative on Q with domain 
consisting of those f G Cq{Q) which satisfy condition (1.1b). 

Remark 3.3. Note that condition (1.1a) entails that if iv^ and /3 are defined by (3.5) 
then Wk G [0, 1], = 1, . . . , n, Wk = 1, and (3 > 0. Therefore the choice (3.5) is 
consistent with the conditions on these parameters required by the construction of the 
elastic Walsh process W^. 

Next we compute the resolvent of the elastic Walsh process. As a byproduct 
this will give another proof of theorem 3.2. Moreover, it will provide us with an 
explicit formula for the scattering matrix in this case. In contrast to the calculations 
in [15, Chapter 2.3], [18, Chapter 6.2] for the classical case with Q = M+, we do 
not use the first passage time formula (1.26), but instead we use formula (1.37). This 
simplifies the computation considerably. 
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Let / G Co{Q), A > 0, and ^ G ^. In the present context formula (1.37) reads 

Rim) = Rj:f io - exio RViv), 

where B}" is the resolvent of the Walsh process without killing, and e\ is defined 
in (1.15). The Laplace transform of the density of C/3,o under Py is readily computed: 

Lemma 3.4. For all A, /3 > 0, 

-AC/3,0 ) ^ 



Proof. As remarked before, we may consider U" to be equal to the local time at the 
origin of the Brownian motion B underlying the construction of W, and therefore the 
analogous statement is true for the right continuous pseudo-inverse K'^ of L^. As 
above let denote the right continuous pseudo-inverse of (cf. 1.2). Then for 

s e M+, 

= e-^^ 

where we used lemma 1.2. Hence 

Jo 

which proves the lemma. □ 
With lemma 3.4 we obtain the following formula 

(3.6) Rifio = Rj:m) - exio Rj:f{v). 



Note that R^f is in the domain of the generator of the Walsh process, and therefore 
satisfies the boundary condition (2.2): 



Y,MRxf)\vk) = Q. 



k=l 



On the other hand, we obviously have e\{vk) = — v 2A for all k G {1,2,..., n}. 
Thus with Ylk=i Wfc = 1 we find. 



J2wk{Rlf)\vk) = P-^j=R^f{v), 
k=l P + VZA 



while equation (3.6) yields for ^ = v 



\/2A 

' /3 + V2A ' 
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The last two equations show that for all / G Co(^), A > 0, we have 
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k=l 

Since for every A > 0, i?^ maps Co{Q) onto the domain of the generator of W^, we 
have another proof of theorem 3.2. 

Upon insertion of the expressions for the resolvent kernels of the Walsh process, 
equations (2.3), and (2.4), with the same notation as in lemma 2.3 we immediately 
obtain the following result: 

Lemma 3.5. For X > 0, 7] & G the resolvent kernel of the elastic Walsh process 
is given by 



k,m=l 
n 



P + V2X 



(3.7b) 

with the scattering matrix S' 
(3.7c) 5|^(A) = 2 



1 

(ri), 



k,m=l 



2A 



P + V2\ 



Wm-Skm, X> 0, k,m e {1,2,. .. ,n}. 



Remark 3.6. Note that in contrast to the case of the Walsh process, this time the scat- 
tering matrix is not constant with respect to A > 0. Also, when /3 = 0, formula (2.4c) 
is recovered, as it should be. In analogy with the discussion in remark 2.4, the bound- 
ary conditions for the elastic Walsh process is given by the matrices 
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such that 



S^{\) = Sa-,b4E = -2A) 

= -{A" + V2XB^)-^ {A' 



2XB^). 



Observe that for k, m G {1, ... ,n} the matrix element ^^^(A) of the scattering 
matrix is obtained from the resolvent kernel as 



5^m(A) = V2A lim {rl{^,v)-rx{^,v)), 
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= ^ (E S^^iX) + l) for all ke{l,2,..., n}, 



where the limit on the right hand side is taken in such a way that ^, ij converge to v 
along the edges l/-, Im respectively. S^^{\) in turn fixes the data Wra and (5, e.g., via 
the behavior for large A, that is from the behavior at "large energies" 

Wm = \ {hm + lip 'S'fcm(-^')) > ^r all G {1, 2, . . . , n}, 

Z A'too 

and 

= V2A( ^^- + ^^7J-f";^^^^ -l), forallA,andalU,mG{l,2,...,n}. 

Alternatively, the data can be obtained from the small A behavior, that is the threshold 
behavior, of the scattering matrix, since from 

-f = lim — = {St^iX) + 6km) for all G {1, 2, . . . , n}, 

we obtain 

2V2A 
and therefore 

li.nH.A-/^(5|JA) + M faall*,*'e{l,2....,n}. 

lim«„A-i/2(E,„SJ,.„(A) + l) 

Furthermore we remark that in the context of quantum mechanics in the self-adjoint 
case Wk = l/n, k = 1, . . . , n, the boundary conditions of the elastic Walsh process 
are interpreted as the presence of a (5-potential of strength P at the vertex. 

In order to compute expressions for the transition kernel of the elastic Walsh pro- 
cess, we use the following two inverse Laplace transforms which follow from formu- 
lae (5.3.4) and (5.6.12) in [6] (cf. also appendix C in [20]) (A > 0, i > 0, x > 0): 

(3.9) ^—e-^^- ^ g^t,x)-^ef'-+^''/'evk(^+/3^[l). 

P + V2X ' 2 V 2/ 

Then the inverse Laplace transform of the scattering matrix is given by the follow- 
ing measures on M+: 

s%^{dt) = {2wm - Skm) eo{dt) - 2wmf3 dt 

(3.10) ^ 

+ Wml5^eP""^ erfc(/3^^)dt, 

with k,mE: {1,2,..., n}. Moreover, for t > 0, x > 0, let us introduce 

(3.11) g0,o{t,x) = git,x) - \ e^^^+^^^Z^ erfc(-^ + /J^l). 
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Lemma 3.7. For t > 0, ^, rj E Q, the transition density of the elastic Walsh process 
is given by 

n 

(3.12) p'{t,^,v)=p''{t,^,v)+ E h,{O^^m90fl{t^d,{^,v))liJv), 

k,m=l 

and alternatively by 

n pi 

(3.13) k,m=l 

d{v,r,)))it-s))lUr,), 

where * denotes convolution. 

4. The Walsh Process with a Sticky Vertex 

In this section we construct Brownian motions on Q with a = in the boundary 
condition (1.1). 

Consider the Walsh process on ^ from section 2 together with a right continu- 
ous, complete filtration J^'^, relative to which it is strongly Markovian. Furthermore, 
we denote its local time at the vertex v by L"" (cf. section 3). 

Again we follow closely the recipe given by Ito and McKean in [14] (cf. also [18, 
Section 6.2]) for the case of a Brownian motion on the half line. For 7 > introduce 
a new time scale r by 

(4.1) T'^ : t ^ t + jL'l' , t>0. 

Since L"^ is non-decreasing, is strictly increasing. Moreover, we have t~^(0) = 
and lim^^+oo T~^{t) = +00, which implies that r exists, and is strictly increasing 
from M+ onto M+, too. As is shown in [18, p. 160], the additivity of L"' entails the 
additivity of r on its own time scale, i.e.: 

Lemma 4.1. For all s, t>0, a.s. the following formula holds true 

(4.2) r{s + t) = T{s)+T{t)o9^^,y 

It is easily checked that for every t > 0, T(t) is an J^^'-stopping time, and since r 
is increasing, we obtain the subfiltration = {F^, t > 0) of J-"^ defined by = 
F^^^y t e R+. Moreover, we set F^ = cf{Ff , t G M+) and F^ = aiFf, t G M+), 
and find F^ C F^. Standard calculations show that the completeness and the right 
continuity of F'^ entail the same properties for F^. (For details of the argument in 
the case where Q = M+ we refer the interested reader to section 3 of [20].) 

Define a stochastic process on Q, called Walsh process with sticky vertex, by 

(4.3) Wi' = Writy teR+. 

Observe that when W is away from the vertex, is constant, and therefore in this 
case r^^ grows with rate 1. On the other hand, when W is at the vertex, r^^ grows 
faster than with rate 1, and therefore r increases slower than the deterministic time 



22 



V. KOSTRYKIN, J. POTTHOFF, AND R. SCHRADER 



scale t ^ t. Thus "experiences a slow down in time" until W has left the vertex. 
In this heuristic sense the vertex is "sticky" for W^^, because it spends more time there 
than 

Note that because has continuous paths, and therefore also r are pathwise 

continuous. Consequently, has continuous sample paths. Since W has contin- 
uous paths, it is a measurable process, and hence for every t > 0, W^^f^ is ^^{ty 
measurable, that is, is J"*-adapted. Set 6*1 = 6*^(4). With the additivity (4.2) of r 
we immediately find 

(4.4) W! o 01 = W!+t, s, t G M+. 

Thus 6^ = {9f, t e M_|_) is a family of shift operators for W^. 

Next we show the strong Markov property of relative to following the 
argument sketched briefly in section 6.2 of [18] for the case Q = M+. First we prove 
the simple Markov property of W'"^ with respect to J^*. To this end, let s, t > 0, 
4 G t?, and C € 8(0). Then we get with (4.4) 

= ^«(^r(s)Oe.(,)GC|j-,-(,)) 

= Pwi{w:eC), 

where we used the strong Markov property of W with respect to T^. As a next 
step we prove that has the strong Markov property for its hitting time of the 
vertex. By construction, and W have the same paths up to the hitting time of 
the vertex, and in particular is also the hitting time of the vertex by W, that is, 
= H^. Moreover, since L^{H^a) = 0, we get that t~^{Hv) = = T{Hy), as 
well as 9'{Hy) = 9{Hy). Assume now that t > 0, ^ G ^, and C G ^(e?). Then on 
{Hy < +00} we can compute with the strong Markov property of W as follows 

PdwUH. e c I ^]^J = p^{wi oef,^€C\ j-]^J 
= PdWrit)oeH^ec\j^^J 

= P.(T^.(i) gC) 
= P,{WfeC). 

It is readily checked that J^f^^ C J^fj^ , and therefore we get in particular the strong 
Markov property of with respect to = Hy in the form 

(4.5) P^ {WUhs G I J = n {Wi G C) . 

Finally, with the strong Markov property of the standard one-dimensional Brownian 
motions on every edge and the strong Markov property (4.5) just proved we can apply 
the arguments of the proof of theorem 1.4.3 to conclude that is a Feller process. 
Hence it is strongly Markovian relative to the filttation J^^. 
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By construction, is up to time equivalent to a standard one-dimensional 
Brownian motion, and it has continuous sample paths. Hence, altogether we have 
shown that is a Brownian motion on G in the sense of definition 1.3.1. 

Now we want to compute the generator of W^, and first we argue that v is not 
a trap for W^. To this end, we may consider W as constructed from a standard 
Brownian motion B as described in section 2. Let Z denote the zero set of B. Given 
s > we can choose to > s in the complement Z*^ of Z. Consider t = T~^{to), 
i.e., t = to + jLf^y Obviously t > s, and r(t) G Z"^. Therefore -B^(t) ^ 0, and 
consequently = Vl/V(t) / ^• 

Theorem 4.2. Consider the boundary condition (1.1) with a = 0, c G (0, 1), and 
b G [0, 1]". Set 

(4.6) Wk = - ,k = l,...,n, 7=- , 

1 — c 1 — c 

and let be the sticky Walsh process as constructed above with these parameters. 
Then the generator ofW^ is 1/2 times the second derivative on Q with domain 
consisting of those f G Cq{Q) which satisfy condition (1.1b). 

Before we prove theorem 4.2 we first prepare two preliminary results. Let e > 0, 
and let ^ denote the hitting time of the complement of the open ball B^{v) with 
radius e and center v by W^. Recall that if^^ denotes the corresponding first hitting 
time for the Walsh process W. 

Lemma 4.3. P^-a.s., the formula 

(4-7) i?,% = i/- +7^^,^, 

holds true. 

Proof. Let W , and therefore also W^, start in the vertex v. Since and W have 
continuous paths with infinite lifetime we have for all 7 > 

= inf > 0, d{v,W^^i)) = e], 

and in particular for 7 = 0, 

i^^, = mi{t > 0, d{v, Wt) = e}. 
Moreover, as argued above, both infima are a.s. finite. Set 

_ ttW I ^, rw 

Then T{a) = -ff^^, and therefore 

d{v,W^)=d{v,W^^^)) 
= d{v,WHwJ 
= e. 
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Consequently we get H^^^ < a. To derive the converse inequality we remark that 

e = d{v,W!jsJ 

which impUes 

Since r is strictly increasing tiiis entails 

and the proof is finished. □ 
Corollary 4.4. For every 7 > 0, 

(4.8) E^{HQ=e^ + je 

holds. 

Proof. By construction, the paths of W starting in v hit the complement of B^{v) 
exactly when the underlying standard Brownian motion B (cf. section 2) starting at 
the origin hits one of the points ibe on the real line. Thus under P^, L'^{H^^^) has 
the same law as L^{H^_^ under Pq. Lemma 1.3 states that under Pq this random 
variable is exponentially distributed with mean e. Then equation (4.8) follows directly 
from lemmas 4.3, and 2.1. □ 

Given these results, we come to the 

Proof of theorem 4.2. Let Wk, k = 1, . . . ,n, and 7 be defined as in (4.6), and note 
that due to the condition (1.1a) on bk, k = 1, . . . , n, and c, we have Wk G [0, 1], 
k = 1, . . . , n, J2k '^k = 1> as well as 7 > 0. Hence we can construct the associated 
sticky Walsh process as above. 

Let denote the generator of with domain V{A^). Then we have for / G 
P(^*), A^f{v) = l/2f"{v) (cf. lemma L5.2). On the other hand, we can compute 
A^f{v) viaDynkin's formula as follows 

EJf{W%H^^,)))-fiv) 
A'f{v) = lim ^ ^ 



= lim 



6^0 E^{HQ 



€4,0 + 76 ' 

where we used corollary 4.4. Since the directional derivatives of / at 

/>fc)= lim kG{l,2,...,n}, 
a^v,a&ii, a[a,v) 

exist (cf. section 1.5), we obviously get the boundary condition 

(4.9) lf"iv)=^j2''knvk) 

^ k=l 
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as a necessary condition. Finally, inserting of the values (4.6) of the parameters w^, 
k = 1, . . . ,n, and 7 into equation (4.9) we complete the proof of theorem 4.2. □ 

Next we shall compute the resolvent of the Walsh process with sticky vertex. 
Similarly to the alternative proof of theorem 3.2 for the elastic Walsh process, as a 
byproduct we obtain an alternative proof of theorem 4.2. We begin with the following 

Lemma 4.5. Let X> 0, f £ Cq{Q). Then 

(4.10) liRVYiv) = ^\ (2A(e^,/) - V2Xf{v)) 



'2\ + 7A 
holds, where 

(4.11) e'^{0=wkex{0, ^ e k, k = 1, . . . ,n, 

and ex is defined in equation (1.15). 

Proof. Let be the generator of PF* on Co (^). From the identity = Ai?^-id, 

and the definition of r we get 

\ {Rif)"{v) = \E, [j^ e-'\f{Wn - f{v)) dt) 

= XE, e-^(^+^^") ifiWs) - f{v)) {ds + jdL^)) 
= \E, e-Mm^n ^f^Ws) - f{v)) ds) . 

In the last equahty we used the fact that only grows when W is at the vertex v. 
By construction of the Walsh process W we have 

E,(e->^^^^f{Ws)-f{v))) 

n 

= Y,^kEo{e-^^''f{fk{\BS - A(0))) 
k=l 

= 2E^'^y^ I e-'^y{Mx)-MO))^==e-(^+y^y'^dxdy, 

where we used lemma 1.1. We insert the last expression above, and use formula (1.31). 
This gives 

-| ^ 1 poo 

^ ^ V2A + 7A 7o 

J— (2A(e^,/)-V2A/(.;)). □ 

2A + 7A 



From the identity A^Ry^ = XR\ — id and some simple algebra we get the 
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Corollary 4.6. Let A > 0, and f e Co{Q). Then 

(4.12) Rifiv) = \ (2 (e^, /) + jf{v)) 

v2A + 7A 

holds. 

Since formula (1.27) in corollary 1.5 is valid for the resolvent of every Brownian 
motion on Q, we may use that formula for R^^f, sum it against the weights Wk, k = 
1, . . . , n, and insert the right hand side of equation (4.12). This results in 

n 

k=i V 2A + 7A 

and a comparison with formula (4.10) shows that equation (4.9) holds true for / re- 
placed by i?^/ for arbitrary / G Co{Q). As promised we thus have another proof of 
theorem 4.2. 

With the help of the first passage time formula we can now provide explicit expres- 
sions for the resolvent R^, its kernel and the transition kernel of W". Inserting 
the right hand side of equation (4.12) into the first passage time formula (1.26), we 
immediately obtain for / G Co(^), A > 0, 

(4.13) Rim = R^f{0+ eA(e)(2(e^,/) + 7/(^)), ^eQ, 

V 2A -|- 7A 

where R^ is the Dirichlet resolvent (1.22). Using formula (1.23) for the kernel of R^ 
together with (1.26), and (1.27), we get the following result. 

Corollary 4.7. For ^, tj £ G, X > 0, the resolvent kernel r^, of the Walsh process 
with sticky vertex is given by 

n ^ 

riiC^ drj) = rf (C, rj) drj + V ex,kiO 2ti^m -7= r eA,m(??) drj 

1 ; , V2A + 7A 

(4.14) k,m=l ^ ' 

+ ^^r^ — reA(0 ^v{dv), 

V 2A -I- 7A 

with defined in (1.23), and denotes the Dirac measure in v. Alternatively, is 
given by 

n ^ 

rxiCdrj) = rx{C,v)dv+ Yl ^\k{0 SlmW ^=e\,m{ri) dr] 

(4.15a) k,rn=l 

+ /^^, ^ ^x{0^v{dn), 

V 2A -I- 7A 

where rx is defined in equation (1.24), and 



(4.15b) Si^{\) = 2 Wm - 5km. 

V 2A -I- 7A 
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Remark 4.8. When all Wm, m = 1, . . . , n, are equal to 1/n, the matrix ^^(A) takes 
the form 

*''^^^-' + vi^^" 

which reduces to (2.9) when 7 = 0. S^{X) is unitary for all A < 0. Also the S'*(A) 
for different A all commute. As a consequence ^^(A) has the interpretation of a 
quantum scattering matrix in the sense of [23]. More precisely, S'^{X) stems from 
the Schrodinger operator —A*, where is a self-adjoint Laplace operator on L?'{Q) 
with boundary conditions of the form (2.5) with the choice 

A = -\[S\\,) + 1), 

(4.16) 1 

B = --^{S%\o)-l), 



for any Aq for which \/2\q + 7A0 7^ 0. We emphasize that the Schrodinger operator 
—A* and the generator of the Walsh process are quite different: Not only do 
they act on different Banach spaces, but also the functions in the intersection of their 
domains satisfy different boundary conditions at the vertex v. As matter of fact, the 
integral kernel of the resolvent (— A* + 2A)~^ of the Schrodinger operator —A* is 
given by, see Lemma 4.2 in [24], 

1 " 1 



k,m=l ^ 

that is — up to a factor 2 — by the right hand side of (4.15a) without the last term. 
In more detail and with the definition (2.6) 

Sa,b{E = -2\)=S'{\) 

holds for all A > 0. As a function of k {k'^ = E), is meromorphic in the complex 
A;-plane with a pole on the positive imaginary axis at k^ = 2i/^. This corresponds 
to a negative eigenvalue = —4/7^ of —A*. The corresponding (normalized) 
eigenfunction ip^ — physically speaking a bound state — is given as 



2 y n 



So quantum mechanically the vertex v acts like an attractive potential. We view this 
as a quantum analogue of the stickyness of the vertex v. 

This analogy can be elaborated a bit further by inspecting the associated quantum 
mechanical time delay matrix (see, e.g., [1,4, 16,25-27,32]) 

which in the present context gives 

^^^^ " M4 + Zv) 
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So T{k) has zero as an (n — l)-fold eigenvalue plus the non-degenerate eigenvalue 

-27 
A;(4 + A;V)' 

which for 7 > is the signal for a strict quantum delay. Observe that for k — > +00, 
that is for large energies, the time delay experienced by the quantum particle tends to 
zero, while for /c — 0, i.e., for low energies, the delay becomes arbitrarily large. From 
the physical point of view, both effects are clearly to be expected. For comparison 
and in contrast to the present stochastic context, in quantum mechanics 7 < is also 
allowed for a meaningful Schrodinger operator and an associated scattering matrix. 

Define for x > 0, 7, t > 0, 

. X 1 /2x 2t\ / X V2i\ 

(4.17) go,^[t,x) = - expl h erfcl ^= H . 

It is not hard to check that 

(4.18) \imgo,j{t,x) = g{t,x) = e'^'/^*. 
74-0 v27rt 

Moreover, from [6, eq. 5.6.16] (cf. also appendix C in [20]) the Laplace transform is 

(4.19) jC,go,^{-,x){X) = -=^ -e-^^ x>0. 

V2A + 7A 

Observe that in agreement with (4.18) 

ji:g{.,x){X) = ^e-^'^ 



holds. Now we can readily compute the inverse Laplace transform of formulae (4. 14), 
(4.15), and obtain the following result. 

Corollary 4.9. For t > 0, ^, rj E G the transition kernel of the Walsh process with 
sticky vertex is given by 

n 

(4.20) + Yl \iO^Wmgo,^{t,dy{C,il))li^{r])dri 

k,m=l 

+ 19o,j{t,d{^,v)) e^{dr]) 

where is defined in equation (1.21), or alternatively by 

p''{t,^,dvi) = p{t,^,r])dr] 

n 

+ Yl \iO{^'(^m9o,j{t,dy{^,r)))dr) 

(4.21) k,m=l 

- Skm9{t,dy((,,r])) dr]^ h^iv) 
+ 19o,-y{t,d{(,v)) e^idr]), 
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andp{t, a, h) is given in formula (1.17). 

We close this section with some remarks concerning the local time of at the 
vertex v, which also serve to prepare the construction of the most general Brownian 
motion on the single vertex graph Q in the next section. 

Let us define 

(4.22) Lt = L^it), t>0, 

where — as before — L'^ denotes the local time of the Walsh process at the vertex, 
having (cf. section 3) the same normalization as the local time of a standard one- 
dimensional Brownian motion (cf. (1.8)). By construction, is pathwise continuous 
and non-decreasing. It is adapted to T^, and a straightforward calculation based on 
the additivity of L"' and formula (4.2) shows the (pathwise) additivity property 

(4.23) Ll^t = LI + LI o 91, s,t>0. 

Thus L^ is a PCHAF of {W^,T^). Furthermore, t > is a point of increase for L^ 
if only if r(t) is a point of increase for L"^, which only is the case if is at the 
vertex, i.e., if Wf is at the vertex. Thus, it follows that L^ is a local time at the vertex 
for W^. In order to completely identify it, it remains to compute its normalization, 
and it is not very hard to compute its a-potential (the interested reader can find the 
details for the case Q = M+ in [20]): 

(4.24) eJ [ e""* dL^) = ^-^d{i,v) ^ a>0,^eg. 

Vo ^ V2a + 'ja 

5. The General Brownian Motion on a Single Vertex Graph 

Finally, in this subsection we construct a Brownian motion by killing the Walsh 
process with sticky vertex of section 4 in a similar way as in the construction of the 
elastic Walsh process (cf. section 3). realizes the boundary condition (1.1) in its 
most general form. 

Consider the sticky Walsh process W'^ with stickiness parameter 7 > 0, right 
continuous and complete filtration J^^, and local time L'^ at the vertex. We argued 
in section 4 that L^ is a PCHAF for {W^,T^), and therefore we can apply the 
method of killing described in section 1.5: We bring in the additional probability 
space (M_|_,;S(M+),P;3) where is the exponential law of rate (3 > 0, and the 
canonical coordinate random variable S. Then we take the family of product spaces 
{n,A, {P^, ^ e g)) of {n,A, {P^, ^ e g)) and (M+,H(M+),P^). Define the ran- 
dom time 

(5.1) = iiif{t > 0, L| > S}. 

Then by the arguments given in section 1.5, the stochastic process defined by 
Wf = Wf for t G [0, C/3,7)' and = A for i > (f^^^, is again a Brownian motion 
on g in the sense of definition 1.3.1. 

Denote by the right continuous pseudo-inverse of L*. Since is continuous 
(cf. equation (4.22)), we get L|^s = r for all r G M+. Recall that the right continuous 
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pseudo-inverse of the local time U" of the Walsh process was denoted by K'^ . Then 
we have the following 

Lemma 5.1. For all 7 > 0, the following relation holds true: 

(5.2) K'^ = + 7r, r G M+. 
Proof. For 7, r G M+ define the random subset 

Mr) = {t>0, LI > r} 

of M+. Since is pathwise increasing, this set is a random interval with endpoints 
and +00. The relation = r implies that 

J-y(r) = {K^,+oo). 
In particular, we have Jo = {K^, +00). Now 

t G J^{r) <^ LI = > r r(t) G Jo(r). 

In other words, J'^'{r) = T^^(Jo(r)), and therefore = T~^{Kf) holds. From 
the definition of r^^ (see equation (4.1)), and the relation = r we obtain for- 
mula (5.2) ' □ 

In the proof of lemma 3.4 the Laplace transform of the density of Kf, r > 0, under 
Py has been determined as A exp(— \/2A?')- Hence we have 

As a consequence we find the 

Corollary 5.2. For r >0, has the density 

(5.3) Py(K' G = , ^ g-rV2(/-7r) ^ > 

Furthermore, the Laplace transform of the density ofK^ under P„ is given by 

(5.4) £;^(e-^^^^) =e-(^+^^)^ A > 0. 

Remark 5.3. One can use lemma C.l in [20] to check that the right hand side of 
equation (5.3) is indeed the inverse Laplace transform of the right hand side of for- 
mula (5.4). 

Observe that (p^j = Kg and = Kg. Thus we obtain the 
Corollary 5.4. For all /3 > 0, 7 > 0, the following equation holds true 
(5-5) C/3,7 = C/3,0 + 7'S'- 

As before, denotes the expectation with respect to P^, ^ G Q. 
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□ 



Corollary 5.5. For a// /3 > 0, 7 > 0, A > 0, the following formula holds true 

(5.6a) 4(e-^^'^'^) =MA), 

with 

(5.6b) p{\) = . 

Proof. With corollary 5.2 and C/3,7 = Kg we obtain 

j'OO 

Jo 

Denote by the resolvent of W^. With lemma 1.6 we immediately find the 
Corollary 5.6. For all f € Co{Q), X > 0, ^ G Q the following formula holds true: 

(5.7) RlfiO = RlfiO - PpiX) exiO Rlf{v). 

Now it is easy to verify that for appropriately chosen parameters /3, 7, w^, k = 
1, . . . , n, the Brownian motion Wg realizes the boundary condition (1.1b). 

Theorem 5.7. Consider the boundary condition (1.1), and assume that b is not the 
null vector Set r = a + c G (0, 1), and 

(5.8) Wk = - ,k = l,...,n, 13 =- , 7 



17-- J 7 ■ - 7 I- 1 7 I -I 

— r 1 — r 1 — r 

Let be the Brownian motion as constructed above with these parameters. Then the 

generator ofW^ is 1/2 times the Laplacean on Q with domain T>{^AP^ consisting 
of those f G Cq{Q) which satisfy condition (1.1b). 

Proof. As in the previous cases it is readily seen that the definition (5.8) of the pa- 
rameters 7, /?, Wk, k = 1, . . . ,n, is consistent with the conditions used in the above 
construction of W^. 

Let be the generator of with domain 2?(A^). Since is a Brownian 
motion on Q in the sense of definition 1.3.1, it follows from lemma 1.5.2 that V{A3) C 
CiiQ), and that for all / G V{A9), A9f{0 = 1/2 /"(O, C & G- Let h G Co{g), 

follows 



A > 0. Then RVti G V{A^), and therefore we may compute with equation (5.7) as 



I {Rih)"{v) = I {Rih)"{v)-f5p{X)2X{Rih){v) 



= J2wk {Rih)'{vk) -Pp{X)jX{Rih){v), 



k=l 
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where we used the fact that, since Rj^h is in the domain of the generator of W^, it 
satisfies the boundary condition (4.9). We rewrite this equation in the following way: 

n 

I {Rlh)"{v) = Y^wu {Rlh)'{vu) + pV2Xp{X){Rlh){v) 
(5.9) ^ fe=i 

-/3/9(A)(^ + 7A) {Rlh){v). 

Now we differentiate equation (5.7) at ^ G Z^, A; = 1, . . . , n, let ^ tend to v along any 
edge Ik, and sum the resulting equation against the weights w^, k = 1, . . . ,n. Then 
we get the following formula 



(5.10) {RihYivk) = Y.Wk {Rih)'ivk) + f3V2Xp{\){Rih){v), 

k=l k=l 

where we used Ylk '"^fc ~ ^- ^ other hand, for ^ = v, equation (5.7) gives 

(5.11) {Rlh)iv) = pi\){V2X + 7A) {Rih) (v) . 

A comparison of equations (5.10), (5.11) with (5.9) shows that we have proved the 
following formula 

(5.12) I {Rih)"{v) = Y,Wk {Rlh)'{vk) - P{Rih){v). 

k=l 

With the values (5.8) for /3, 7, and Wk, k = 1, . . . , n, it is obvious that / = 
satisfies equation (1.1b). Since i?^ is surjective from Cq{Q) onto the domain of the 
generator of W^, the proof of the theorem is finished. □ 

Let A > 0, / G Co(^). Insertion of the right hand side of formula (4.13) for R^ 
into equation (5.7) gives us after some simple algebra the following expression for 

(5.13) RlfiO = R^fiO + P(A) e^iO (2(e^, /) + ^f{v)) , ^eQ, 

where R^ is the Dirichlet resolvent, e\ is defined in equation (1.15), in equa- 
tion (4.11), and p{\) is as in formula (5.6b). From equation (5.13) we can read off 
the following result: 

Corollary 5.8. For ^, f] E Q, X > 0, the resolvent kernel of the general Brownian 
motion on Q is given by 

n 

riii, dv) = rf (C, v)d'n+ V ex,k{0 2w^m p(A) ex,m{v) drj 

+ 1pW exiO ^vidrj), 
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with as in formula (1.23), and p is defined in equation (5.6b). Alternatively, 
can be written in the following form 

n ^ 

. ''a(^' di) = rx{^, ri)dr}+ ^ ex,k{0 Sl^{X) -i= ex,m{r]) drj 
(5- 15a) V2A 

+ 7P(A) exiO ev{dri), 

where rx is defined in equation (1.24), and 

(5.15b) 5^^(A) = 2 V2Ap(A) - Skm- 

In order to invert the Laplace transforms in equations (5.14), (5.15), we define for 
^, 7 > 0, the following function g/s^^ on (0, +oo) x M+: 

with (t, x) G (0, +oo) X M_|_. The heat kernel 5^,^ is discussed in more detail in 
appendix C of [20]. In particular, it is outlined there that the limits of gp^j as /3 4 
0, and 7^0, yield the kernels 5^3,0 (equation (3.11)) and ^0,7 (equation (4.17)), 
respectively. Moreover, it is proved there that the Laplace transform of gp^j{-,x), 
a; > 0, is given by 

(5.17) /9(A) e-^^, A>0, 
where p is defined in (5.6b). Hence we get the 

Corollary 5.9. For C> ^ ^ 0, t > 0, the transition kernel of the general Brownian 
motion on Q is given by 

p^{t,^,dr]) =p^{t,C,v)dv 

n 

(5.18) + hdO^Wmg0,j{t,d^{^,r,))liMdri 

k,m=l 

+ 19lS,-f{t,d{^,v)) ev{dr]), 
which alternatively can be written as 
p^{t,^,dr)) =p{t,^,r])dr] 

n 

m 5/3,7 

(5.19) k,m=l 

~ Skmg{t,dyiS,,T]))^ li^ir])dr] 
+ 19/S,t{t,d{^,v)) ey{dri). 
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